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Abstract 


Non-classical stress concentration behavior in a stretched circular hyperelastic sheet (outer radius 
b= 10 in., thickness ¢ = 0.0625 in.) containing a central hole (radius a = 0.5 in.) was analyzed. The 
hyperelastic sheet was subjected to different levels of remote radial stretchings. Nastran large-strain large- 
deformation analysis and the Blatz-Ko large deformation theory were used to calculate the equal-biaxial 
stress concentration factors K. The results show that the values of K calculated from the Blatz-Ko theory 
and Nastran are extremely close. Unlike the classical linear elasticity theory, which gives the constant 
K =2 for the equal-biaxial stress field, the hyperelastic K values were found to increase with increased 
stretching and can exceed the value K = 6 at a remote radial extension ratio of 2.35. The present K-values 
compare fairly well with the K-values obtained by previous works. The effect of the hole-size on K- 
values was investigated. The values of K start to decrease from a hole radius a = 0.125 in. down to K = 1 
(no stress concentration) as a shrinks to a = 0 in. (no hole). Also, the newly introduced stretch and strain 
magnification factors {K,,K,} are also material- and deformation-dependent, and can increase from 
linear levels of {1.0, 4.0} and reaching {3.07, 4.61}, respectively at a remote radial extension ratio of 
2.39; 


Nomenclature 

a hole radius, in. 
b remote radius, in. 

A, 0 0 
Cy =| 0 A 0 |, Cauchy-Green deformation tensor in principal axes 

0 0 A 

WwW — 
Cc’ = ar” Mooney-Rivlin constant 

1 

WwW a 

CS = ap Mooney-Rivlin constant ( J, -dependent) 
2 
d mathematical symbol for differentiation 
E Young’s modulus, Ib/in? 
re =TrC, =C, =A; +A; + Ai, first invariant of C, 
1 
I, = »CiCz = s(n — Cia) Ap AS + AAS +AZA;, second invariant of C,, 
i#k 
f, =DetC,, =A; A5A; , third invariant of C, 
I =[,=C,=A, +A +A, first invariant 
J at (G5 : + : + : modified second invariant 
2 =e Ee ii 427 42° 42? 
: Ay as AS 
Deformed volume 
J =,/1, =./DetC, =1,A,A, = ———_—_———— modified third invariant 
: : x" Undeformed volume 

i,j = 1,2,3, indices 
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Oo 


max “max 


= am , stress concentration factor 


0... co 
A ie de 
= 7 , Stretch magnification factor 
E ; : ‘ 
=— , strain magnification factor 
E 


radial coordinate, in. 
thickness, in. 


unmanned aerial vehicle 
strain energy density function, in-Ib/in* 


ow 
= va = C,, first strain energy gradients (material constant), Ib/in? 
1 
OW 2 
= a7. = C,J,, second strain energy gradient (material constant), Ib/in 
2 
ow 


= ——, third strain energy gradient, Ib/in” 
3 


=A —-1, strain, in/in 

equal-biaxial strain, in/in 

maximum tangential strain at hole boundary, in/in 

uniaxial strain, in/in 

tangential strain, in/in 

biaxial strain at remote boundary, in/in 

extension ratios in 1, 2, 3-directions 

= 1+ €, extension ratio in loading direction 

extension ratio in lateral directions for uniaxial tension 
extension ratio in thickness direction for biaxial tension 

= dp(r)/dr, radial extension ratio 

= p(r)/r tangential extension ratio 

=/,(b) = (dp/dr),_,, remote (r = b) radial extension ratio 
=1,(a)=(p/r),_,, hole-boundary (r = a) tangential extension ratio 
shear modulus at small strain (A, > 1), Ib/in? 

= In Aree 


= , modified Poisson’s ratio for large deformations 


Ind 


deformed length of r, in. 


uniaxial engineering stress (Biot stress) in reference to an undeformed cross section, 


Ib/in? 


strip-biaxial engineering stress in reference to an undeformed cross section, Ib/in 


equal-biaxial engineering stress in reference to an undeformed cross section, lb/in 


engineering stress in a radial direction, 1b/in 


engineering stress in a tangential direction, Ib/in? 


O.. =0,(b), remote (r = b) radial engineering stress, Ib/in? 

Ownax =0,(da), hole boundary (r = a) tangential engineering stress, Ib/in? 

O, radial true stress (Cauchy stress) in reference to deformed cross section, Ib/in? 

O ini =(0,,,A)/J,, uniaxial true stress (Cauchy stress) in reference to deformed cross section, 
Ib/in? 

O, tangential true stress (Cauchy stress) in reference to deformed cross section, b/in’ 


Introduction 


Hyperelastic (rubber-like) materials have wide applications (for example, automobile tires, rubber 
hoses, shock absorbers, inflatable unmanned aircraft, hovercraft skirts, et cetera). In the major aerospace 
applications, rubber materials are traditionally used as fuels in solid rocket motors (for example, Space 
Shuttle solid rocket boosters, surface-to-air, and air-to-air missiles). Rubber burns very slowly, however, 
when mixed with oxidizer particles (for example, ammonium perchlorate oxidizer particles) to form 
composite solid propellants, the solid rocket propellants can burn very rapidly to generate tremendous 
thrust. For example, each Space Shuttle solid rocket booster could generate approximately 3,300,000 Ib 
thrust at lift-off. Because of particle inclusions, the composite solid propellant has stress concentration 
problems. For example, the star tips of the star-grain (fig. 1), when ignited, could be the potential crack 
initiation sites because of tensile stress concentration induced by the chamber pressure. 

Because of increasing aerospace applications of hyperelastic materials [For example, transition flaps, 
launch vehicle engine vibration damping, structural flutter suppressions, inflatable wing skins of an 
unmanned aerial vehicle (UAV), et cetera.], there is a need to fully understand fundamental mechanical 
behavior of such materials. The hyperelastic membrane has no compressive strength. When applied to 
inflatable aircraft skins, or transition flaps for noise reduction (ref. 1), the hyperelastic skins must be pre- 
stretched to maintain shapes. When a sharp particle strikes and penetrates a stretched hyperelastic sheet, a 
small hole generated could cause a very high stress concentration around the hole boundary. An excessive 
high stress concentration can cause instant tearing failure. A pressurized balloon is a good example of a 
highly stretched hyperelastic skin. With a simple contact with a sharp object, the balloon can burst 
instantly because the pinhole created induces a high stress concentration at the pinhole boundary causing 
the hole boundary material to reach failure stress. 

In general, uniaxial stress-strain data are used as a basis to describe the mechanical behavior of 
hyperelastic materials. However, in aerospace applications, hyperelastic skins are usually under biaxial 
stress fields. Therefore, there is a need to fully understand the biaxial stress-strain behavior and non- 
classical stress concentration behavior of the stretched hyperelastic materials. 

In the classical linear elasticity theory, the stress concentration factor K around a hole is K = 3 under 
the uniaxial stress field, and K = 2 under the equal-biaxial stress field. Those K values are independent of 
material properties and independent of stretching (ref. 2). However, for large deformations of a 
hyperelastic sheet with a hole, the value of K can be dependent of material properties and hole-size, and 
can increase with increased stretching. 

This technical report will analyze the non-classical stress concentration problem in a radially 
stretched circular hyperelastic sheet containing a central circular hole made of widely-used Mooney 
material (refs. 3, 4), and will theoretically confirm that the equal-biaxial stress concentration factor K is 
material, hole-size, and stretch dependent. 

Nastran (ref. 5) large-strain large deformation analysis (using fine mesh) and the Blatz-Ko large 
deformation theory (refs. 6, 7) were used to calculate the hyperelastic equal-biaxial stress concentration 
factors K. The results show that the values of K calculated from Nastran and from the Blatz-Ko large 
deformation theory (refs. 6, 7) are extremely close, and compare fairly well with the K-values obtained in 
other analyses (refs. 8-11). 


Unlike linear elasticity theory, which gives the constant K = 2 for the equal-biaxial stress field, the 
equal-biaxial stress concentration factor K for hyperelastic materials was found to be material-dependent, 
hole-size-dependent, and deformation-dependent, and can increase with increased stretching, reaching the 
values several times that of K = 2 (linear case). 

Also, in the report, the stretch- and strain-magnification factors for hyperelastic materials are 
introduced, and the material- and deformation-dependent behaviors are discussed. 


Material Considered 


This section describes the historical Mooney material, which was used in the current biaxial stress 
concentration analysis, and presents some comments about using the reduced plot to obtain hyperelastic 
material constants. 


Mooney Material 


For the current rudimentary stress-concentration analysis, the historical Mooney material (W, = 
18.35 Ib/in’, W, = 1.835lb/in’) (refs. 3, 4) was chosen. The two material constants {W,,W,} are the two 
types of strain energy gradients defined as W,=0W/0J,=0W/dl,=C, ; J,=1,=4; +25 +A; and 
W, =0W/dJ, =(OW/O1,)J,=CyJ,3 J, = Ap tVAp+i/AZ; 1, =A; +A +A; 

J,= JI, = DetC, =1,1,,. 

Figure 2 shows the uniaxial stress-strain curve of the Mooney material, and figure 3 shows the 
rey) 
18.35 lb/in’, W, = 1.835lb/in’) Ib/in’ through linear curve fitting. Because the Mooney material was 


reduced plot, against 1//, using the data of figure 2 for obtaining the material constants ( W, = 


widely used by earlier researchers (refs. 8—11) in the similar biaxial stress concentration studies, by using 
the same Mooney material in the current biaxial stress concentration analysis, the results can be easily 
compared with the earlier results (refs. 8-11). Therefore, one can examine the consistency of the biaxial 
stress concentration behavior of the same hyperelastic material. 


Comments on the Reduced Plot 


uni 


PVA against 1/A shown in figure 3, 


It is important to mention that in the reduced data plotting, 


the value of Paar at the limit A=1(o 


=0) is not zero nor indeterminate (0/0), but has a finite 


uni 


value equal to shear modulus w [=(2W,+2W,)] of incompressible materials (v = 0.5). Namely, by 


expanding 


V-VWA V-1/Aa 


will become shear modulus fu of incompressible materials (v = 0.5) according to the following expanded 


; oA 
Suni h_ in the neighborhood of A +1 (that is, € 0), the limit value um 
Jz! 


equation (1): 
Frith Belle) atl Je(I+e) 3 
ON EU) = 8 Ee) 2 EG 
Rah. 2 Geer ate Che -daer) Be 


In general, the reduced plot (similar to figure 3) for any hyperelastic material is nonlinear, but 
contains a large linear region. For the linear curve fitting to determine the material constants 
{2W,,2W,} , one can use the value of mu calculated from equation (1) to estimate the location of the 
vertical intercept (2W,+2W,) at 1/A=1 of the curve-fitting straight line (similar to figure 3). 


As an example, using the value of elastic modulus E =117.92 lb/in’, determined from the initial 


Oni 
slope of the uniaxial stress-strain curve of figure 1, the limit value (f= | (= LL) for the 
- del 
; : ; EB L792 4 
incompressible Mooney material can be calculated as: H= = —— = 39.31 lb/in’, 
DEV) le 3 


which is extremely close to the value of (2W,+2W,)= 40.37 lb/in’ (vertical intercept at 1/A =1) 
determined from linear curve fitting shown in figure 2. 


Descriptions of the Problem 


For the stress concentration analysis, the problem considered is a circular hyperelastic sheet 
containing a central circular hole subjected to remote radial stretching (fig. 4). Figure 4 shows the 
geometry of the hyperelastic sheet with a boundary radius of b = 10 in., thickness ¢ = 0.0625 in., and hole 
radius a = 0.5 in., and is made of Mooney material (W, = 18.35, W, = 1.835) lb/in’. Figure 5 shows the 
Nastran (ref. 5) finite-element quarter model for the case b = 10 in. The case of b = 5 in. was first 
considered, however, an accurate remote equal-biaxial stress field could not be obtained because b = 5 in. 
was too small. Therefore, b had to be increased to b = 10 in. in the current analysis. At b = 10 in. the 
remote boundary will be practically under an equal-biaxial stress field. The hyperelastic sheet was 
subjected to different levels of remote radial stretching (A,, = 1.25, 1.50, 1.70, 1.98, 2.10, 2.35) (used by 
earlier authors, refs. 7-10). Nastran large-strain large-deformation analysis was then used to study the 
effect of stretching on the stress concentration factor K. Also, the Blatz-Ko large deformation theory 
(ref. 6, 7) was used to calculate the stress concentration factors K and compared with the corresponding 
Nastran-calculated K values. 


Deformation Equations 


If radial length r of the hyperelastic circular sheet with a hole is radially deformed to p(r), the radial 
and tangential extension ratios {A, Ao} can be expressed respectively with equation (2) 
(refs. 2, 3): 


2,()= PO 5 Airy=2O 2) 


At the deformed hole boundary, the tangential extension ratio i, (a) can be written as equation (3): 


(3) 


At the deformed remote boundary, the radial and tangential extension ratios {A. (b),A, (b)} become 
identical if b is sufficiently large (b = 10 in.). Namely, one can establish equation (4): 


A... =40)-(2 ) = Z ia =A,(b) (4) 


The remote radial extension ratio A,, is to be specified, but the tangential extension ratio A at 


deformed hole-boundary must be calculated from equation (3) by inputting the experimentally measured 
deformed hole radius p(a) . 


For a special case when the hole shrinks to zero, and the specimen degenerates into a continuous 
sheet, the radial and tangential extension ratios become identical as shown in equation (5): 


G2 BO ©) 


=A,(r) (5) 
Thus, the strain field becomes the equal-biaxial strain field everywhere in the continuous specimen. 


Large Deformation Theory 


This section presents the general form of the Blatz-Ko stress-difference equation (refs. 6, 7), and its 
special forms written respectively for the uniaxial, strip-biaxial, and equal-biaxial stress fields for 
applications to the calculations of the biaxial stress concentration factor. 


Stress-Difference Equations 


By taking the difference between stresses in 7 and j directions, one can eliminate the terms containing 
W,(=0W/dJ,) in which J,= (i= {DetC,, =A,A,A,. The resulting Blatz-Ko stress difference 
equation has the following form shown in equation (6) (refs. 6, 7): 


0,A,-90 A, 2W, 7 - 
“-e =2W, + eA ; (i,j=1,2,3) (i,j not summed) (6) 


In equation (6), {0,,0,}and{A,,/,} are respectively the engineering stresses (Biot stresses) and 
extension ratios in (i,j) (=1,2,3) directions, and {W,(=0W/dJ,);W,(=0W/dJ,)} are the material 
constants. For incompressible materials (J, = de =1,1,A,=1), the quantities {(0,A,),(o jAp sin 


equation (6) will become the true stresses (Cauchy stresses) respectively in (i,j) (=1,2,3) directions. 


Equation (6) is applicable to both compressible (for example, foamed rubbers) and incompressible 
hyperelastic materials. The stress-difference equations (6) can now be written for three major stress fields 
(uniaxial tension, strip-biaxial tension, and equal-biaxial tension) for incompressible materials (J, = 1) 


as presented below. 


I. Uniaxial Tension 


6, =6,=02A4,=A, A, =A, =4, 


For uniaxial tension (O, =O lat 


), the stress difference equation 


uni ? 
2 
lat 


(6) for incompressible materials (J, =/,A,A, = AA,,, =1) takes on the form of equation (7) (refs. 6, 7): 


O nih 2W. 
uni = 2W. + 2, 
LA Oe a 7) 
By plottin Ow against 1/A for the uniaxial case, one can curve-fit the data points with a linear 
yp 8 42 8 


trend line, and determine 2W, from the vertical intercept at 1/A =O and (2W, +2W,) from the vertical 
intercept at 1/A = | (fig. 3). 


2. Strip-Biaxial Tension 

For strip-biaxial tension (0,=0,,;,0,40,0,=0; A,=A,A,=1,4,=A,), the stress 
difference equation (6) for incompressible materials (J;=A,A,A,=AA,,=1) takes on the form of 
equation (8) (refs. 6, 7): 


Ch 
age ee (8) 
By plotting o,_,,A against (A° os 1/ A’), one can curve-fit the data points with a linear trend line and 
determine the summation of (2W,+2W,) from the slope of the plot, but 2W, and 2W, cannot be 
determined separately from the plot. Keep in mind that in the biaxial balloon test, the stress field near the 


ballooned specimen boundary will be under the strip-biaxial stress field, and can be described by 
equation (8). 
3. Equal-Biaxial Tension 

For the equal-biaxial tension (0, =O, =0,_,,, 0;=0; 4,=4, =A, A, =1,,), the stress difference 


equation (6) for incompressible materials (J,=A,A,A, =A°A,, =1) takes on the form of equation (9) 
(refs. 6, 7): 


O._pi A 


Equation (9) can be rewritten in a more convenient form as equation (10): 


O,,A 2W, 
= +2W. 10 
a ey aaa (10) 
By plotting poe against 1/ A’ , one can curve-fit the data points with a linear trend line, and 


determine 2W, from the vertical intercept at 1/A7=0; obtain (2W,+2W,)= from the vertical 
2 1 =~ 


intercept at 1/ A? =1, and obtain 2W, from the slope of the curve-fitting straight line, or calculate from 
2W,=(W,+2W,)-— 2W, . 
nKH— 
Determined Determined 
For a circular continuous specimen under radial biaxial tension, the entire specimen will be under the 
equal-biaxial stress field, which can be described by equation (9) or (10). However, in the biaxial balloon 
test, the equal-biaxial stress field occurs only in the ballooned specimen center. 


Stress Difference Equations for the Current Problem 


For a radially stretched circular hyperelastic sheet with a central circular hole, the remote boundary 
will be under the equal-biaxial stress field, but the traction free hole boundary will be under the uniaxial 
stress field in a tangential direction (fig. 6). Therefore, the uniaxial stress-difference equation (7) and the 
equal-biaxial stress-diffence equation (9) can be written for the hole boundary and remote boundary for 
stress calculations. 


At the Hole Boundary 


At the hole-boundary, the uniaxial stress difference equation (7) can be written in the following form 
for the Mooney material (W) = 18.35, W2 = 1.835) Ib/in? (refs 2, 3) as shown in equation (11): 


1 
0,4 = 0,(a)A,(a) a Ciodicg = 2 W, + Me Uae = — 
True stress Aroax Arca 


U3=1) (11) 


1.835 1 
= 2] 18.354+— a 
Ags | oe i) 


In equation (11), O,,[=O,(a)] and A,,,.[=,(a)] are respectively the tangential engineering stress 


and the tangential extension ratio at the hole boundary. To calculate the true stress (0,,,An.,) from 


max” “max 


equation (11), the unknown A 


max 


can be determined from the equal-biaxial tensile tests by taking the ratio 
of the deformed hole-radius to the undeformed hole-radius according to equation (3). However, for the 


current analysis, values of A 


max Were estimated from the Nastran large-strain large-deformation analysis. 


At the Remote Boundary 


At the remote boundary, the equal-biaxial stress difference equation (9) can be written for the 
Mooney material (W, = 18.35, W2 = 1.835) Ib/in’ (refs 2, 3) in the following form as shown in equation 
(12): 


6,4 = 6, (b)A,(b)= 0A. = 2(W, + Wo) 23 7 a 


4 
True stress A. 
Us=)) (12) 
= 2(18.35+1.89522)[2 7: 


In equation (12), o,[=o0,(b)] and A,[=A,(b)] are respectively the radial engineering stress and the 
radial extension ratio at the remote boundary. Because A, is specified, the true stress o,A,, at the remote 
boundary can be calculated from equation (12). 


Stress Concentration Factor 


The equal-biaxial stress concentration factor K is defined in equation (13): 


Hole boundary true tangential stress _6,(a)_0,(a)A,(a)_ 0,4 


max “max 


K = = = 
Remote boundary true radial stress 06(b) ob)A(b) oA, us) 


For classical linear elasticity, equation (13) will give a deformation-independent value of K = 2. 
Note that for the current radial loading, if b 210 in., the remote boundary will be under the equal-biaxial 
stress field, and therefore, the radial and tangential stresses and extension ratios become identical; 


namely, 0,(b)=0,(b) =0., and A,(b) = A,(b) = A., . In the current analysis, the outer radius b = 10 in. 


was used for the circular finite-element model (fig. 3) for achieving the equal-biaxial stress field. By 
plotting K as a function of A,,, one can geometrically observe how the theoretically calculated and the 


Nastran calculated values of K increase with increasing /,,. Note that, to obtain the theoretical values of 


K, one has to calculate the true stresses {(0 ),(o,A,,)} respectively from the large deformation 


max Mesa 


constitutive equations {(11), (12)}. 


Stretch Magnification Factor 


For large deformations, measuring the extension ratio is much easier than measuring stresses because, 
no accurate large-deformation stress gages exist. Therefore, the stretch magnification factor K, can be 
introduced to study how the stretch magnification factor K, changes with stretching. Thus, one can avoid 
using large deformation constitutive equations. In the experiments or finite-element analysis, the remote 
extension ratio Ais specified. However, /,,,. must be measured or calculated using finite-element 
analysis. Once ines is determined, then one can calculate the stretch magnification factor K, from 
equation (14): 


_ Hole-boundary tangential extension ratio _A,(a) _ A 


max 


K = = 
i‘ Remote radial extension ratio Afb) dz uy) 


‘co 


Then one can plot the curves {K, against 2} and study how the value of K, changes with increasing 


A,,. When the strains approach zero, K, will approach the limit case of K, =1. 


Strain Magnification Factor 


It is also a good idea to use the strain magnification factor K, and explore how the strain 
magnification factor K, changes with stretching. When Arnaz iS Obtained from measurements (or finite- 


max 


element analysis), then one can calculate the strain magnification factor K, from equation (15): 


max — Emax 


. Remote radial strain A(b)-1 €,(b) A,-1 «€, KP?) 


KR Hole boundary tangential strain A,(a)—-1_ €,(a)_A,,,, —1 


Then, one can plot the curve of K,as a function of A,,for comparisons with the limit case of linear 
theory. 


Limit Cases of the Strain Magnification Factor Based on the Large Deformation Theory 


For small deformations, the uniaxial strain € ) and the equal-biaxial strain €.(=€,,.) can 


(= E e—bi 
be calculated respectively from the large deformation stress equations (7) and (9) by expanding in the 
neighborhood of A—1. Keep in mind that equations (7) and (9) are written for incompressible 


materials. 


max uni 


1. Limit Uniaxial Stress at the Hole Boundary 
The large-deformation uniaxial stress equation (7) can be rearranged and expanded in terms of small 
uniaxial strain € (<<1) as equation (16): 


unt 


2W. 1 2W. 1 
Oo. =0,.=| 2W,+— || A-— |—— ] 2W, + 2 l+e,.)-—-————— 
max uni ( 1 yi it a AO 1 1 l+e it ae (l+e,,,)° : 


uni 


=(2W, + 2W,(1-€,,, +...) [(1+ €,,,)-(—26,,; +.) (16) 
= (2W, + 2W, )3€, 4, = 3UE sn 
Tee 


From equation (16), uniaxial strain can be expressed in terms of uniaxial stress as equation (17): 


Eni a — (17) 


2. Limit Equal-Biaxial Stress at the Remote Boundary 


The equal biaxial stress equation (9) can be rewritten and expanded in terms of small strain 
E,_,;(<< 1) as equation (18): 
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O,.9, =(2W, + 2WoA')( 2 -z:] 


= [2W, ze 2W, (I+ €.5:) I +E, 4; eget OG Ee pi ae ] 


(18) 
=[(QW, + 2W,)+2W,(2é, ,,+.-JIA+é,_,, -1+5€,_,, +...) 
= (2W, + 2W, )6€,_,, = OME,_5, 
=H 
Equation (18) can be rewritten in a final form as equation (19): 
€ zs 0 4 19 
e-bi 6u ( ) 


3. The Resulting Limit Case of the Strain Magnification Factor 


By substituting equations (17) and (19) into equation (15), the value of K, for the linear case can be 
calculated as equation (20): 


K, = E max = E wni = O uni 6u = 2 Oni = 2 O max = 4 
E.. Ebi 3u O 2; Oj 0. (20) 
Euni 1/€,_p; =K=2 


Equation (20) shows that under the equal-biaxial stress field, for the linear-case of incompressible 
materials, the strain magnification factor K, =4 is twice the linear-case of the stress concentration factor 


K=2. 
Limit Case of the Strain Magnification Factor Based on the Linear Theory 


For small strain, the limit case of strain magnification factor K, given by equation (15) can also be 
calculated using the linear theory as shown in equation (21): 


Uniaxial strain at hole boundary 


1 
bec = nS OUT SoG. HO) 
€ Ex €,(d) be ={6,(6)-Vi5,(b)+ 01} 0 .(b)—vo,(b) (l—v)o..(b) 


6, (b)=09(b) 


(21) 


Equal-biaxial strain at remote boundary 
Equation (21) shows that the value of K, is a function of Poisson’s ratio v. For incompressible materials 


(for example, hyperelastic materials), the Poisson’s ratio is practically v = 0.5, then equation (21) takes on 
the form of equation (22): 
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= 0,(a) _—— 1 o@)_ 1 [On oy 
(l-v)o,(b) *y2os —(0.50,(b) O5\ 0, 


v—dependent =K=2 


(22) 


Thus, for incompressible materials (v = 0.5), both the linear theory and the limit case of the large 
deformation theory give the identical equal-biaxial strain magnification factor of K, = 4 [eqs. (20), (22)]. 


Numerical Results 


For the calculations of {K,K,,K,} [eqs. (13), (14), (15)], values of {A,,,A,,,.} are needed. The 
remote extension ratio A,, is specified and is known. However, the unknown A, [= p(a)/a] must be 
determined from experiment. Namely, by simply measuring the deformed hole radius p(a) and taking 
the ratio p(a)/a to obtain A 


simpler than measuring the large deformation stress o,,,. at the small hole boundary because no 


max according to equation (3). The geometrical measurement of A,,,. is much 
miniature strain gages can operate under very large deformations. In the present report, values of y ie 


were estimated by using Nastran large-strain large-deformation analysis. 


Nastran Data 


The major data generated from Nastran large-strain large-deformation analysis are listed in table 1. 
The Nastran-generated K listed in table 1 will be plotted in the subsequent section (Plots of K Curves), 
after the corresponding theoretical values of K are calculated for comparisons. 


Table 1. Nastran generated data for calculating the equal-biaxial stress concentration factor K, and the 
stretch and strain magnification factors {K nN ,K “ay : Mooney material (W, = 18.35, 


W>= 1.835) lb/in?; (a = 0.5, b = 10, t= 0.0625) in. 


X., o.,, Ib/in” oe On, > Ib/in® K K. K, 
1.25 48.26 2.06 141.71 2.94 1.65 4.24 
1.50 90.78 3.09 338.69 3.73 2.06 4.18 
1.70 128.61 3.92 547.36 4.26 2.31 4.17 
1.98 192.41 5.18 956.96 4.97 2.62 4.28 
2.10 224.62 5.79 1192.26 5.31 2.76 4.35 
2.35 302.75 7.22 1850.38 6.11 3.07 4.61 


Average = 4.3048 


Figure 7 shows plots of the Nastran-generated tangential extension ratios Xp (r) as functions of radial 
distance r under different levels of the remote radial extension ratios A,,. Note that for each level of A, 
the corresponding Xp (r) increases toward the hole and reached the maximum value Xe (a) =A,,, at the 
hole-boundary. 

Figure 8 shows the Nastran-generated hole-boundary tangential extension ratio Ass [= A,(a)] plotted 


as a function of the remote radial extension ratio A,,[=A,(b)]. Note that the A, 


max CUrve is practically a 


straight line up to A, =1.75, and then slightly bend upward as A, further increases. At the remote 
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stretching of A, =2.35, the hole-boundary extension ratio A, can reach as large as A, 


(fig. 7). 


Theoretical Predictions of the Stress Concentration Factor 


7.22 


aX 


In the calculations of theoretical stress concentration factors K, the known remote extension ratios J, 


and the Nastran-estimated hole boundary extension ratio 
(Cauchy) radial stresses oA, and hole-boundary true (Cauchy) tangential stresses o 


max 


results of calculations are listed in table 2 with Nastran data (table 1) included for comparison. 


Amax Were used to calculate remote true 


A... The 


max 


Table 2. Calculations of stress concentration factors K using equations {(11), (12)} compared with 


Nastran K; Mooney material (W) = 18.35, W2 = 1.835) Ib/in’; (a=0.5,b 


10, ¢ 


0.0625) in. 


Equal-biaxial for remote boundary 


Uniaxial for hole boundary 


1 1.835 1 
©, A, = 2(18.35 + 1.835 x A2)} AZ-— O iniAmox = 2] 18.354 Rox 
a. Ana Anax 
True stress A, True stress 


B., Amax 


|< pene nen ne nnn nn nae Equal-biaxial at remote boundary--------------- >|< ween enn n nnn n enna anne Uniaxial at hole boundary----- 
2 2 
X.. XQ: Ag B co 0. Ag Dvds yee Le Bi O nax Aves 
[= A, (b)| (Equal-biaxial) [= As (a)| 


(Uniaxial) ( _ Ono Ainax 
~ OLA, 


1.01 1.0201 40.7438 0.0591 2.4088 1.04 1.0816 40.2288 0.1201 
1.25 1.5625 42.4344 1.1529 48.9226 2.06 4.2436 38.4816 3.7582 
1.50 2.2500 44.9575 2.0525 92.2753 3.09 9.5481 37.8877 9.2245 349.4951 
1.70 2.8900 47.3063 2.7703 131.0526 3.92 15.3664 37.6362 15.1113 
1.98 3.9204 51.0879 3.8553 196.9591 5.18 26.7289 37.4099 26.5355 
2.10 4.4100 52.8847 4.3586 230.5033 5.78 33.4084 37.3349 33.2354 
2.35 5.5225 56.9676 5.4897 312.7350 7.22 52.1284 37.2083 51.9899 


4.8315 2.0058 
Nastran = 2.0000 
144.5088 2.9538 
Nastran = 2.9364 
3.7875 
Nastran = 3.7309 
568.7319 4.3397 
Nastran = 4.2560 
992.6904 5.0401 
Nastran = 4.9735 
1,240.8403 5.3832 
Nastran = 5.3079 
1,934.4558 6.1856 
Nastran = 6.1119 


Plots of the Stress Concentration Factor 


Figure 9 shows the plots of theoretical K listed in table 2 and Nastran-generated K of table 1. Note 
that, the two K-curves calculated from both the Blatz-Ko large deformation theory and from the Nastran 
large-strain large-deformation analysis are extremely close. Note also that, unlike the classical linear 
theory which gives material and deformation independent K = 2, the value of K for hyperelastic material 
increases almost linearly (with slight waving) with increased stretching, and can exceed K = 6 at the 


remote extension ratio of A_ = 2.35. 
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Figure 10 compares the K-curves based on current analysis and K-curves obtained by earlier 
researchers. Note that Yang’s K-curve [theoretically calculated using (W,/W,)=0.1 (ref. 8)], which 


only ends at (K =5.4,A,, =2.1), is very close to the K curves calculated from the Blatz-Ko theory and 


the Nastran large-strain and large-deformation analysis. The Segal-Klosner’s K-curve [experimentally 
obtained using a square sheet (6.5 in.x 6.5 in.) of natural rubber (ref. 9)] is a straight line only up to K = 4 
at the remote extension ratio A, =1.62. The last data point falls on the vicinity of the K-curves 


calculated from Nastran and the Blatz-Ko theory (fig. 10). However, Oden’s K-curve (refs. 10, 11) is 
somewhat off and gives higher K values. Keep in mind that Oden’s K-curve was finite-element- 
generated, and the discrepancy could be attributed to the use of a much coarser finite-element model with 
a smaller remote radius of b = 5 in. (ref. 11), and the use of a slightly lower value of (W, /W,)=0.08. In 


the current analysis, it was found that the value of K was extremely sensitive to the finite-element mesh 
sizes and the dimension of remote radius b. 


Plots of the Stretch and Strain Magnification Factors 


Figure 11 shows the stretch magnification factor K, (data from table 1) plotted as a function of the 
remote extension ratio 7, . Note that, the value of K, starts from K, =1 and increases convex-upwardly 
with increasing A,,, and then reached the value of K, = 3.07 at A, = 2.35. 

Figure 12 shows the strain magnification factor K, (data from table 1) plotted as a function of the 
remote extension ratioA,. Note that, the value of K, starts from the linear value of K, = 4, and 


increases nonlinearly with increasing A,,, and reached the value of K, = 4.61 at A, = 2.35. 


Numerical Results of Hole-Size Effect on the Stress Concentration Factor 


After the above stress concentration studies, additional studies were conducted to investigate the hole- 
size effect on the stress concentration factor K in the equal-biaxial strain field. The hole-sizes considered 
are a = (0.5, 0.25, 0.125, 0.05) in. 


Equal-Biaxial Stress Concentration Factor 


The theoretically-calculated K-values and corresponding Nastran generated K-values are listed in 
table 3. 
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Table 3. Theoretical equal-biaxial stress concentration factors K [eqs. (11), (12)] compared with Nastran- 
generated K ; Mooney material (W; = 18.35, W2 = 1.835) Ib/in?; [a = (0.5, 0.25, 0.125, 0.05), b = 10, 


t = 0.0625] in. 


True stress 


OA. = 2(18.35 + 1.835 x a 
Cetin 


1 
at) oA 


uni” “max 


True stress 


= a{18.38 +— 


1.835 1 
An 
Armax I aa Ainax 


A. oA. Ans O inax max K K AK 
(Equal biaxial) (Uniaxial) (Ospax max )/(O-nAce) (Nastran) (% difference) 
a=0.5 in. 

1.01 2.3910 1.04 4.8299 2.0200 2.0000 1.00 
1.25 48.9226 2.06 144.5088 2.9538 2.9364 0.59 
1.50 92.2753 3.09 349.4951 3.7875 3.7309 1.52 
1.70 131.0526 3.91 568.7319 4.3397 4.2560 1.97 
1.98 196.9591 5.17 992.6904 5.0401 4.9735 1.34 
2.10 230.5033 5.78 1240.8403 5.3832 5.3079 1.42 
2.35 312.7350 7.22 1934.4558 6.1856 6.1119 1.21 

a= 0.25 in. 
1.01 2.3910 1.04 4.8299 2.0200 2.0000 1.00 
1.25 48.9226 2.04 141.3455 2.8892 2.9600 2.39 
1.50 92.2753 3.11 354.1997 3.8386 3.7800 1.55 
1.70 131.0526 3.90 562.8685 4.2950 4.3100 0.35 
1.98 196.9591 5.18 996.5380 5.0596 5.0600 0.01 
2.10 230.5033 5.81 1253.7462 5.4392 5.4200 0.35 
2.35 312.7350 7.28 1966.6484 6.2885 6.2500 0.62 

a=0.125 in. 
1.01 2.3910 1.04 4.8299 2.0200 1.9900 1.51 
1.25 48.9226 2.01 136.4813 2.7897 2.8700 2.80 
1.50 92.2753 3.03 335.5472 3.6364 3.6500 0.37 
1.70 131.0526 3.86 549.7804 4.1952 4.1500 1.09 
1.98 196.9591 5.13 977.3639 4.9622 4.8800 1.68 
2.10 230.5033 5.66 1189.8801 5.1621 5.2200 1.11 
2.35 312.7350 TAT 1907.8306 6.1005 6.0100 1.51 

a = 0.05 in. 
1.01 2.3910 1.03 3.6247 1.5160 1.6000 5.25 
1.25 48.9226 1.83 108.4702 2.2172 2.1500 3.13 
1.50 92.2753 2.65 253.0796 2.7427 2.6700 212 
1.70 131.0526 3.33 407.8317 3.1120 3.0200 3.05 
1.98 196.9591 4.35 701.7895 3.5631 3.5000 1.80 
2.10 230.5033 4.84 869.7430 3.7733 3.7300 1.16 
2.35 312.7350 6.06 1363.8403 4.3610 4.2600 2.37 
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Plots of Hole-Size Effect Data 


Figure 13 shows the plots of A,,,, as functions of A,, for different hole sizes. The A,,,, -curves for 


the hole-sizes a = (0.5, 0.25, 0.125) in. almost form a single curve, indicating that under those hole-sizes, 
there is practically no hole-size effect on the stress concentration factor K. The A,,,, -curve for the hole- 
size a=0.05 in. is somewhat off from the group of A,,,, -curves, and has a lower slope, showing the 
effect of the hole-size on K. 

Figure 14 shows the stress concentration factors K plotted as functions of A,, for different hole sizes. 


The K curves for hole-sizes a = (0.5, 0.25, 0.125) in. are very close, indicating that under those hole- 
sizes, there is practically no hole-size effect on K. The K-curve for the hole-size a = 0.05 in. is 
considerably off from the group K-curves, and has lower slope, indicating the existence of hole-size effect 
on K. 

Figure 15 shows the hole-boundary tangential extension ratio A,,,, plotted as functions of hole- 


radius a for different levels of remote extension ratios 2,,. As shown in table 2, A.,,, is the key control 


factor affecting the stress concentration factor K, and therefore, it is important to first examine the 
behavior of the A.,,, -curves (fig. 15) before the K-curves are presented. As shown in figure 15, the A,,,, 


-curves for different levels of A,, are practically horizontal in the hole-size range of (0.125 < a < 0.5) in., 


indicating that there is no hole-size effect. When the hole-size is decreased from a = 0.125 in., the value 
of Ajnax for each case of A, starts to decrease with a decreasing-hole radius a, and finally down to A,,,, 
= 1 ata =0 in. (no hole). The pictorial shapes of A,,,, suggest that the corresponding K- curves will have 
very similar geometrical shapes. 

Figure 16 shows the equal-biaxial stress concentration factors K plotted as functions of hole-radius a 
for different levels of remote extension ratios A,,. Note that, the geometrical shapes of the K-curves are 


very similar to those of the corresponding A,,,,-curves shown in figure 12. The K-curves in the range of 
(0.125 < a < 0.5) in. for all levels of A, are practically horizontal, indicating that there is no hole-size 
effect on K. When the hole-size is decreased from a = 0.125 in., the value of K for each case of A, starts 


to decrease with decreasing hole radius a and finally down to K = 1 (no stress concentration) at a = 0 in. 
(no hole). Thus, for the current special case of Mooney material (W; = 18.35, W2 = 1.835) lb/in’, the 
hole-size effect on K will appear only in the small region (0 < a < 0.125) in. 


Concluding Remarks 


1. For the linear elasticity theory, the equal-biaxial stress concentration factor K is K = 2, and is 
independent of strains and material properties. 


2. The strain magnification factor K,for linear case is K,=4(=2K) for incompressible 


materials. 


3. For large deformations, the equal-biaxial stress concentration factors K calculated from the Blatz- 
Ko large-deformation theory and from the Nastran large-strain large deformation analysis using 
Mooney material (W, = 18.35,W, = 1.835) Ib/in? are extremely close, both of which compare 
fairly well with the K-values obtained by earlier authors. 
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For the hyperelastic Mooney material considered, the equal-biaxial stress concentration factor K 
is a strong function of stretching and material property. At the remote extension ratio of 
A, = 2.35, the value of K can exceed K = 6. 


For the hyperelastic Mooney material considered, the equal-biaxial stretch magnification factor 
K,, is a function of stretching, and can increases from K, =1 (limit case) up to K, =3.07 at 


the remote extension ratio of A, = 2.35. 


For the hyperelastic Mooney material considered, the equal-biaxial strain magnification factor 


K, is a weak function of stretching, and can increase from K,, = 4 (linear case) up to K, = 4.61 


at the remote extension ratio of A,, = 2.35. 


For the hyperelastic Mooney material considered, the hole-size effect on the equal-biaxial stress 
concentration factor K starts to appear when the hole radius a is less than a = 0.125 in. Namely, 
at the hole radius a = 0.125 in. the value of K (associated with any remote extension ratio /,, ) 
starts to decrease gradually and finally down to the lowest limit value of K = 1 (no stress 
concentration) as the hole-radius a shrinks to zero (a = 0 in.). 
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Figure 1. Tensile stress critical regions in a solid rocket star grain due to internal pressure and inclusion 
effect. 
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Figure 2. Uniaxial stress-strain curve of the Mooney material (W, = 18.35, W, = 1.835) Ib/in’ used in the 
analysis. 


2W,+2W, 
= 40.37 psi 
2W, = 
36.70 psi 
-e- Mooney material 
-e Curve fitting 
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
LA 
470049 


Figure 3. Reduced plot of the data of figure 1 for determining Mooney material constants 
(W, = 18.35, W, = 1.835) lb/in’. 
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t = 0.0625 in. 
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Figure 4. Hyperelastic circular sheet containing a central circular hole subjected to remote radial 
stretching (a = 0.5, b = 10, t= 0.0625) in.; Mooney material (W; = 18.35, Wz = 1.835) Ib/in’. 
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Figure 5. Finite-element quarter model for hyperelastic circular sheet containing a central circular hole 
(a= 0.5, b = 10, t= 0.0625) in.; Mooney material (W; = 18.35, W2= 1.835) Ib/in’. 
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Figure 6. In a radially stretched circular hyperelastic sheet with a central circular hole, the outer boundary 
will be under equal-biaxial tension, while the traction free hole boundary will be under uniaxial tension in 
a tangential direction. 
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Figure 7. Radial variations of tangential extension ratios Xo (=p / r) for different remote radial extension 


ratios A,[=(dp/dr),_,]; Mooney material (W, = 18.35, W, = 1.835) Ib/in’; (a = 0.5, b = 10, 
t= 0.0625) in. 
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Figure 8. Plot of a Nastran-generated hole-boundary tangential extension ratio A,,,. as a function of the 


X 


remote radial extension ratio A,,; Mooney material (W, = 18.35, W, = 1.835) Ib/in’; (a = 0.5, 
b = 10, t= 0.0625) in. 
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Figure 9. Comparison of theoretical and Nastran stress concentration curves for Mooney material 
(W, = 18.35, W, = 1.835) lb/in’; (a = 0.5, b = 10, t = 0.0625) in. 
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Figure 10. Comparison of stress concentration curves based on current analysis and earlier authors; Yang 
(ref. 8), Segal-Klosner (ref.9), and Oden (refs. 10, 11) for Mooney material (W,= 


18.35, W, = 1.835) Ib/in’. 
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Figure 11. Nastran-generated stretch magnification curves; Mooney material (W, 
1.835) Ib/in’; (a = 0.5, b = 10, t= 0.0625) in. 
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Figure 12. Nastran-generated strain magnification curves; Mooney material (W, = 18.35, W, = 


1.835) Ib/in’; (a= 0.5, b = 10, t = 0.0625) in. 
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Figure 13. Plots of Nastran-generated hole-boundary tangential extension ratios A,,,, as functions of the 
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remote extension ratio 2,, for different hole sizes; Mooney material (W, = 18.35, W, = 1.835) lb/in’; 
[a = (0.5, 0.25, 0.125, 0.05), b =10, ¢ = 0.0625] in. 
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Figure 14. Plots of Nastran-generated equal-biaxial stress concentration factors K as functions of the 
remote extension ratio A,, for different hole sizes; Mooney material (W, = 18.35, W, = 1.835) Ib/in’; 


[a = (0.5, 0.25, 0.125, 0.05), b = 10, t= 0.0625] in. 
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with hole-radius a under different remote 


extension ratios A,,; Mooney material (W, = 18.35, W, = 1.835) Ib/in’; [a = (0.5, 0.25, 0.125, 0.05), 


b = 10, t= 0.0625] in. 
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Figure 16. Variation of equal-biaxial stress concentration factors K with hole-radius a under different 
remote extension ratios A,,; comparisons of theoretically-generated K-curves and Nastran-generated 


K-curves; Mooney material (W, = 18.35, W, = 1.835) Ib/in? : [a = (0.5, 0.25, 0.125, 0.05), b = 10, 
t = 0.0625] in. 
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